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ABSTRACT 

Context. Dust is commonly present in weakly radio emitting star-forming galaxies and this dust may obscure the signatures of 
accreting black holes in these objects. 

Aims. We aim to uncover weak active galactic nuclei, AGN, in the faint radio source population by means of deep high-resolution 
radio observations. 

Methods. VLBI observations with a world-wide array at unparallelled sensitivity are carried out to assess the nature of the faint radio 
source population in the Hubble Deep Field North and its flanking fields. 

Results. Images of twelve compact, AGN-driven radio sources are presented. These represent roughly one quarter of the detectable 
faint radio source sample. Most, but not all of these low power AGN have X-ray detections. 

Conclusions. The majority of the faint radio source population must be star-forming galaxies. Faint AGN occur in a variety of (distant) 
host galaxies, and these are often accompanied by a dust-obscured starburst. Deep, high-resolution VLBI is a unique, powerful 
technique to assess the occurrence of faint AGN. 

Key words, galaxies: active - galaxies: radio continuum - galaxies: starburst 



1. Introduction 

Deep wide-field radio observations of the Hubble Deep Field 
North (HDF-N) and surrounding Flanking Fields (HFFs) have 
revealed a population of sub-mJy and microJy radio sources, 
of which a substantial fracti on are associated with relatively 
distan t star-forming galaxi es dRichards et al.lll998l iGarrett etall 
|2T3T)11 IMuxIqw et alj|2005h . Optical studies of these fields indi- 
cate that a significant fraction of the hosts of these radio sources 
are faint, remainin g undetected in deep Hubble Space Telescope 
/-band imaging (Richa rds et al.l IT999I) . Some of these objects 
can be seen in the near-IR and have very red colours, lead- 
ing to their interpretati on as distant, dust-enshrouded starbursts 
dDickinson et aljhoOOl) . 

While X-ray and follow-up surveys have established a sub- 
stantial surface den sity of AGN among t h ese and other dis- 
tant galaxies (e .g. lAlexander etall 120031; iBauer et al.l 120041: 
IXue etall 12011 1) , the actual AGN surface (and space) density 
is most likely higher. Many models for the growth of super- 
massive black holes predict that the most rapid growth phase 
will be obscu red by Compton-thick absorption (lFabianl ll999: 
iHopkins etal] f2006). which is for inst ance sup ported b y mid- 
infrar ed observations of Type-2 QSOs (Ma rtmez-Sansigre et al.l 
120051) . The importance of strongly obscured AGN became 
apparent from studies of radio-excess ( with respect to in - 
frared) AGN without X-ray counterparts (iDonlev et al.l 120051) . 
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and from the X-ray properties of seemingly normal or opti 
cally obsc ured star-forming galaxies displaying e xcess infrared 
emission (Daddietal. 2007; Fiore et al. 2009; Murnhv et 
200 91 iPooe et alJ2008b rGeorgaka kiset al.l2010HLuo et al.l201 
Alexande r et al.l 1201 lh . Infrared and submm photometric sur 
veys have proven extremely valuable. The former permit AGN 
identi fication on the basis of the shape of the infr ared SED 
(e.g. lAlonso-Herrero et al.l [20061: IDonlev et aTll2007l). or using 

i; L; 



infrared color-color selec tion (e.g. IStern et al.l l2005: Lacv 
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2004; IDonlev et al.ll2Tn~2n . The latter survey technique finds dis- 
tant starburst galaxies with an occasional embedded AGN (e.g. 
iPope et al.ll2006r. iLutz et"aDl2010t iBiggs et alJl20Tol) . The gen- 
eral reliability of these selection techniques, and the level of 
AGN completeness nevertheless remain debatable. 

Access to the complete infrared-submm spectral energy dis- 
tributions, SEDs, has recently become possible through Herschel 
Space Observatory data, and the SED data permit separation of 
AGN- and starburst-heating in a consistent fashion. On the basis 
of extensive Herschel observations, redshift evolution of star- 
formation in a large sample o f X-ra y selected AGN hosts was 
explained by Mull anev et alj d2012l) as being due to evolution 
of the specific star-formation rate, sSFR, in otherwise normal 
star-forming host galaxies. Employing multi-wavelength SED 
decomposition, a frequent occurence of excess radio emission 
attributed to weak AGN in distant star-formin g galaxies was very 
recently established bv lDel Moro et al.1 (120121) . 

One of the key questions in extragalactic astronomy is the 
true incidence of AGN in distant starburst galaxi es and the 
symbiotic interaction of the two phenomena (e.g. ISilkl [2005: 
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ISchweitzer et al.l 120061) . All above mentioned studies make it 
very clear that the phenomena of star-formation and accretion- 
driven nuclear activity can occur in concert. It is however also 
very clear that dust obscuration plays an important role, pre- 
venting full and complete assessment of the occurrence of AGN 
among mildly or strongly star-forming galaxies using ultraviolet, 
optical, infrared and even X-ray observations. 

Dust obscuration does not affect the radio band, hence cm- 
wavelength radio emission at arcsec scale resolution can be used 
to find (weak) AGN among star-forming galaxies^. The radio- 
infrared ratio can be used to isolate radio-excess o bjects among 
the fa int radio source population, just on itself (e.g. lDonlev et alJ 
120051) , or in combination w ith radio morphology and spectral 
shape dSevmour et al. 2008). Considering the nature of the faint 
radio source p opulation, the power of subarcsec radio imaging 
was shown by Muxlow et al. (2005), and in the case of distant 
infrared galaxies by ICasev et al.l (120091). The present study fol- 
lows up on the work of iMuxlow et al.l (120051) . employing Very 
Long Baseline Interferometry, VLBI, to obtain 1-10 milliarcsec 
scale resolution on the 92 faint radio sources reported by these 
authors in the HDF-N and its flanking fields. Only such deep 
VLBI observations are capable of distinguishing between radio 
emission generated by star-formation processes and AGN activ- 
ity in distant dust-obscured systems. 

Previous VLBI observations of th e HDF-N, carried out with 
the European VLBI Network, EVN dGarrett et al.ll200ll) already 
demonstrated the power of deep, high resolution VLBI imaging 
in discriminating between starburst and AGN activity in dust ob- 
scured systems, yielding important morphological information 
as well as brightness temperatures. This paper presents the cul- 
mination of that study, employing a global array of VLBI tele- 
scopes attaining milliarcsec-scale resolution and micro-Jansky 
sensitivity over a considerably larger field. 

2. Global VLBI observations and data analysis 

Observations of the HDF-N and HFF region were made with 
a global VLBI array on 20-22 February 2004 at 1.4 GHz. A 
total of 36hrs observing time was split into three 12hrmns. The 
global VLBI array consisted of 16 telescopes in Europe and 
the USA, including the 100-m Effelsberg, 100-m Green Bank, 
and 76-m Lovell (Jodrell Bank) telescopes, recording data at 
a bit rate of 128 Mbits/s (2x8 MHz bands) in both right (R)- 
and left (L)-hand circular polarisation. The observations were 
made in phase-referencing mode with a typical cycle time of 
160sec on the secondary phase-calibrator and 330sec on the 
HDF-N. The total on-source integration time on the HDF-N was 
~19 hours. Two phase-calibrators were used: a strong primary 
calibrator J1241+602 (a compact Si.4ghz ~ 455 mJy source 
lying 2° from the HDF-N) and a fainter secondary calibrator 
J1234+619 (S 1,4 ghz ~ 20 mJy lying only 20' from the HDF-N). 
The primary calibrator was typically observed for lOOsec once 
every ~42min; more frequent scans on the nearby, secondary 
calibrator permitted to make the finer phase corrections. 

The data were processed at the European VLBI Network 
(EVN) correlator at the Joint Institute for VLBI in Europe, 
JIVE, in Dwingeloo, the Netherlands. Each 8 MHz band (R&L) 
was correlated in separate passes in order to maximise spectral 
resolution (256 channels/IF) thereby minimising bandwidth 
smearing. The total data set was 675 Gbytes in size making 

1 Radio-loud AGN obviously stand out in radio surveys through their 
radio morphology and luminosity. 



this one of the largest and most complex VLBI data sets ever 
processed. An integration time of 0.25seconds was employed 
to reduce the effects of time smearing. The phase center was 
coincident with a 470yuJy radio source, VLA J123642+621331 
(Richards 2000), located just outside the HDF-N (in an adjacent 
HFF; see Sect. 3. 1.1). The data from each IF were averaged, 
edited, and calibrated with the NRAO AIPS package, and 
gain parameters were subsequently applied to the unaveraged 
original data. The visibility amplitudes were calibrated using 
the system temperatures and gain information provided by each 
telescope. 

In order to image out the entire HDF-N and HFF, the 
wide-field imaging technique was used dGarrett et al.l 1 19991) . A 
large number (92 - see below) of dirty images and dirty beams 
for each 8 MHz band (R&L simultaneously) were generated 
using IMAGR, and co-added; fields with a clear detection 
were subsequently CLEANed using APCLN. The absolute flux 
density scales of the VLBI observations are expected to be 
better than 5%. However, we note that the VLBI flux densities 
have not been corrected for primary beam attenuation. The 
calibrators used have positions that are measured to better than 
one milliarcsec. We therefore estimate the uncertainties in VLBI 
positions to be set by errors introduced by the ionosphere. With 
a separation of 0.5 degrees between our main phase reference 
source and the target fields, we estimate an error of ~1 - 2 
milliarcsec in each coordinate. 



3. Global VLBI detections in the HDF-N and HFF 
region 

The wide-field global VLBI survey covered a total of 201 
arcmin 2 divided into four annular fields with different sensitivi- 
ties and angular resolutions. The observations are both deeper 
and wider than any previous VLBI observations of the field, 
attaining an r.m.s. noise level of 7.3 /Jy/beam with 4 mas 
(milli-arcsec) angular resolution in the central 0-2' part. The 
outer annuli, at 2'-4', 4'-6', and 6'-8' from the phase center, 
have 14-37 yuJy/beam r.m.s. noise level and 19-27 mas angu- 
lar resolution. Radio images were made for all of the 92 ra- 
dio sources which w ere detected by the MERLIN- VLA survey 
dMuxlow et alJl2005l) in these areas. These 92 sources had been 
classified on the basis of their arcsec-scale radio morphologies 
combined with their infrared flux densities, as follows: 18 AGN 
or AGN candidates, 3 starbursts with embedded AGN, 45 star- 
bursts or starburst ca ndidates, and 26 unclassified radio sources 
dMuxlow et alj|2005l) . Taking the variable VLBI sensitivity over 
the whole field into account, 48 out of these 92 radio sources 
are bright enough to be detected by our observations. However, 
VLBI detection critically depends on the angular dimension of 
the radio component(s): extended arcsec-scale emission will be 
completely resolved by the small VLBI beam. 

Our new global VLBI observations clearly reveal 12 com- 
pact radio sources in the HDF-N and HFF above the local 
5cr detection level. This corresponds to a ~25% detection rate. 
Figure 1 presents the radio-optical overlay of the fields, combin- 
ing the arcsec and milliarcsec radio imaging. The global VLBI 
detections are listed in Table 1, with the MERLIN- VLA ra- 
dio structure description and object classification in columns 10 
and 11. The radio structural description is as follows: C: com- 
pact, at 0.2 arcsec resolution; CE: compact with 2-sided ex- 
tended emission; C1E: compact with 1 -sided extended emis- 
sion; FR1: Fanaroff and Riley Class 1 edge-darkened double- 
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lobed radio galaxy. The table also contains ^24-values, defined 
as the logarithm of the 24ji/m/1.4GHz flux density ratio (where 
the integrated arcsec scale radio emission as measured with 
the VLA was used in the computation). Among the VLBI 
sources, three (VLA J123642+621331, VLA J123644+621133, 
and VLA J123646+62 1404) had been det ected in earlier EVN 
1.6 GHz observations (Garr ett et all 12001b . but the other nine 
radio sources are new. All detected objects are described sep- 
arately, in two Subsections below. 

3. 1 . Inner HDF-N and adjacent HFF 

3.1.1. VLA J1 23642+621 331 

VLA J123642+621331 lies just outside the HDF-N, in an ad- 
jacent HFF. It has a steep radio spectru m (a — 0.94, S v oc 
y~ a ); its optic al counterpart has / = 25 dRichards et al.l ll998: 
Mux low et alj|2005l) . HST NICMOS 1.1 yum (F110W), 1.6 fim 
(F160W), and KPNO 4-m /T-band imaging detect a red host 
(^814 — K = 2.0 + 0.2) for the radio source, and spectra obtained 
by the Keck II telescope show a single strong Lya emissio n line 
at z = 4.424 dWaddington et all 1 19991: lDickinsonl[T998h . This 
source is also detected in the ISO 15 urn and Chandra soft X-ray 
bands dAussel et al.lll999HBrandt et al.ll2001l) . IWaddington et all 
(1 19991) interpret the object as a dust-obscured starburst with an 
embedded AGN. 

The deep ME RLIN- VLA observations at 1.4 GHz 
(iMuxlow et al.ll2005l) resolved VLA J123642+621331, showing 
extended emission adjacent to an unresolved core, which corre- 
sponds to the core detected by the EVN 1.6 GHz observations 
dGarrett et al.l 12001). The ratio of its radio and FIR lu minosity 
indica tes that the object has a modest radio excess (Garrett 
2002), providing further evidence for an embedded AGN: its 
1724-value is -0.35. 

The deep, high-resolution global VLBI image removes all 
doubts, revealing a jet-like extension (~ 120 /vJy) emanating 
from a compact A GN core. The integ rated flux density is in 
agreement with the iGarrett et al.l (1200 1 ) measurement. The sep- 
aration between the jet and the AGN core is about 70 pc. This 
may be an example of a high redshift ULIRG in which the high 
star-formati on rate and the efficiency are enhanced by AGN jet 
activity (e.g. lSiikll2005l) . 

3.1.2. VLA J1 23644+621 133 

VLA J123644+621133 is optically identifi ed with an I = 20.9 , 
very red elliptical galaxy at z — 1.050 (ICohen et alj 12000). 
Its large-scale radio structu re and luminosity indicate a clas- 
sical FR-I radio galaxy dRichards et al.1 1 1 9981: IMuxlow et al.l 
12001 . The AGN nature also follows from its nega tive (724 up- 
per lim it. The previous EVN 1 .6 GHz observations dGarrett et alj 
1200 11) revealed a compact core component and an additional 
5cr component located ~60mas south of the core component, 
of which the direction coincid es with the extension of the 
VLA-onl v dRichards et all 19981) and MERLIN- VLA radio mor- 
phology (i Muxlow et alj 120051) . The WSRT and more recent 
dMorrison et al. 20101) VLA flux densities indicate substan- 
tial resolution effects. Our high resolution VLBI observations 
merely detect the core. 



3.1.3. VLA J1 23646+621 404 

VLA J123646+621404 has a slightly inverted radio spectrum 
(a = -0.04) and is associated with a nearly face-on spi- 
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12 36 42.095 42.094 42.093 42.092 42.091 42.090 42.1 
RIGHT ASCENSION (J2000) 
Grey scale flux range= 21.9 100.0 MicroJY/BEAM 
Cont peak flux = 9.7656E-05 J Y- BEAM 
Levs = 7.300E-06 * (-5, -3, 3, 5, 7, 9, 11) 



Fig. 2. VLA J123642+621331. In this and subsequent radio im- 
ages, the clean beam FWHM is specified in the bottom left cor- 
ner. Contours are plotted starting at 3 cr and the relevant noise 
figure can be found below each image. The colour coding is 
specified above each image. 
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Grey scale flux range- 39.0 228.0 MicroJY/BEAM 
Cont peak flux = 2.2806E-04 JY/BEAM 
Levs = 1.300E-05 ' (-5, -3, 3, 5, 7, 9, 11) 
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Fig. 3. VLA J 1 23644+62 1 1 33 



ral g alaxy (/ = 20.7), at a redshift of 0.961 dCowie et al.l 
2004). The MERLIN-VLA image shows a compact component 
and two-sided extended emission, and EVN 1.6 GHz observa- 
tions detect a 4cr radio s ource at this position dGarrett et al.l 
1200 It Muxlow et al. 2005). Its radio variability is evident from 
Richards et al.l d!998[). from the present data, and from recent 
VLA data dMorrison et al.ll2010l) . Those VLA data in fact imply 
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q24 = -0.18. Taken together with the broad emission line char- 
acter (Phill ips et al.lll997b iBrandt et al.lbOQll) . there is no doubt 
about the presence of an AGN. 
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12 36 46.336 46.335 46.334 46.333 46.332 46.331 46.330 
RIGHT ASCENSION (J2000) 
Grey scale flux range- 21.9 220.0 MicroJY BEAM 
Cont peak flux = 2.1780E-04 JY/BEAM 
Levs = 7.200E-06 " (-5, -3, 3, 5, 7, 9, 11, 13, 
15) 



Fig. 4. VLA J123646+621404 



3.1.4. VLA J1 23652+621 444 

VLA J123652+621444 is an inverted spectrum radio source 
(a - -0.12) located in an adjacent HFF and observed to vary 
in intensity on the timescale of months, suggesting the presence 
of an AGN dRichards et al.lll998Ll}^hardst2()0()tlMorrison et all 
120101) . The object displays radio excess, as judged from its 
<724-value. MERLIN-VLA observations reveal a compact core, 
which overlies the nucleus of an / = 18.5 elliptical ga laxy at z — 
0.321, and a one-sided jet -like feature to the east (ICowie et al.l 
l2004t iMuxlow et al.ll2005i) . The high resolution VLBI observa- 
tions merely detect the 0.08 mJy radio core. 

3.2. Outer HFF region 

3.2.1. VLA J1 23608+621 035 

VLA J 123608 +62 1035 is a relatively fiat-spectrum radio source 
(a _ 0.36) hosted by an I = 20.6 galaxy at z = 0.681 
dRichardsll2000T: ICowie et al.ll2004l) . The object is bright in the 
infrared, and falls on the radio-FIR correlation. Its compact ra- 
dio c ore is straddled by tw o-sided radio emission extending over 
^ 1" (IMuxlow et alj2005l) : our VLBI detects a slightly extended 
0. 140mJy radio core. 

3.2.2. VLA J1 23623+621 642 

VLA J123623+621642 has a steep radio spectrum (a = 0.63) 
and its MERLIN-VLA image shows a comp act component 
with weak extende d (sub-arcsec) emission ( Richards! l2000t 
IMuxlow et all [2005). The radio source is identified with an 
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RIGHT ASCENSION (J2000) 
Grey scale flux range= 21 .00 70.00 MicroJY BEAM 
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Levs = 7.000E-06 * f-5, -3, 3, 5, 7, 9, 1 1 , 1 3, 
15) 



Fig. 5. VLA J 1 23652+62 1 444 
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RIGHT ASCENSION (J2000) 
Grey scale flux range- 39.30 78.00 MicroJY BEAM 
Cont peak flux = 7.7795E-05 JY/BEAM 
Levs = 1.310E-05 "(-5, -3, 3, 5, 7, 9, 11, 13, 
15) 



Fig. 6. VLA J 1 23608 +62 1 035 



/ = 23.9 galaxy at z = 1.918; its radio and infrared data indi- 
cate a radio-excess. As inferred from the VLBI observations, a 
substantial fraction of the radio luminosity is contained within 
its milliarcsec-scale radio core. 

3.2.3. VLA J1 23642+621 545 

VLA J123642+621545 is an intermediate radio spectrum source 
(a = 0.50), associated with an / = 20.7 galaxy at a red- 



4 



S. Chi et al.: HDF-N and HFF 



3.2.4. VLA J1 23700+620909 



O 42.76 — 




VLA J 123700+620909 represents similar properties to VLA 
J123642+621331. It has relatively strong radio emission and a 
steep radio spec trum (a = 0.89). Th ere is no optical counter- 
part to R = 27 dRichards et al.ll 19981) . The MERLIN- VLA ob- 
servations detected the source, finding a compact radio com- 
ponent and one-sided emission to the north-east. At the same 
time, deep Spitzer IR AC and MIPS imaging as well as SCUBA 
850 yum observations dPope et al.l20 06) also detected this source, 
GN 16. Given its ERO-character as well as the results of the 
multi-wavelength studies in sub-mm, IR, and radio, this source 
is considered to be a star burst galaxy (SFR ~ 1000 M Q /yr), at 
Zphotom. = L68 dPope et alJl2005ll2006h . The VLBI observations 
find a ~ 0.15 mJy compact radio core in this mild radio-excess, 
q24 = -0.08 object. 



12 36 23.555 23.550 23.545 23.540 

RIGHT ASCENSION (J2000) 
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Cont peak flux = 1.9036E-04 JY/BEAM 
Levs = 1.840E-05 " (-5, -3, 3, 5, 7, 9, 11, 13, 
15) 



Fig. 7. VLA J123623+621642 



shift of 0.857 (Ric hardsl l 2000; Hornschemeier et al. 1I2001I). The 
source was also detected by ISO and Chandra. iMuxlow et alJ 
(200 5]) suggest that the compac t component could be an AGN, 
but lHornschemeier et al.l (l200ll) interpreted the X-ray's as origi- 
nating from a starburst system. The VLBI observations uncover 
a compact component of variable flux density in this infrared- 
luminous object. The object falls on the radio-infrared corre- 
lation using the low-state VLA flux density, but definitely also 
contains an AGN. 
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RIGHT ASCENSION (J2000) 
Grey scale flux range- 48.00 90.00 MicroJY/BEAM 
Conf peak flux = 8.8037E-05 JY/BEAM 
Levs = 1 .600E-05 * (-5, -3, 3, 5, 7, 9, 1 1 , 1 3, 
15) 



Fig. 9. VLA J 123700+620909, GN 16 
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Grey scale flux range- 48.0 290.0 MicroJY/BEAM 
Cont peak flux = 2.8781E-04 JY/BEAM 
Levs = 1.600E-05 " (-5, -3, 3, 5, 7, 9, 11, 13, 
15) 



Fig. 8. VLA J 123642+62 1545 



3.2.5. VLA J 12371 4+620823 

VLA J123714+620 823 is associated with a faint, R = 24.25 
(Cap ak et al.l I2004T) red object. Cowie (2010, priv. communi- 
cation) provided a tentative redshift z = 0.847. With a total 
1 .4 GHz flux density of well ove r a mJy, the object h as a moder- 
ate radio excess, as noted also by Donle vet al.ld2005l) . The VLBI 
radio source is somewhat extended, and the WSRT and VLA 
measurements indicate extended emission on the arcsec and sub- 
arcsec scale. The 21/x Ty 8.4 GHz VLA source J123714+620822 
(Richa rds etall 1998)' is 0.3 arcsec away and its relation to 
J123714+620823 is unclear. 

3.2.6. VLA J1 2371 6+621 51 2 

VLA J123716+6215 12 has a relatively flat radio spectrum (a = 
0.41) dRkhards 2000). The radio structure consists of a compact 
component and one-sided emission extending ~0.4 arcsec to the 
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3.2.7. VLA J1 2371 6+621 733 



62 08 23.35 — 



VLA J123716+621733 is a 0.35mJy radio source with an / = 
22.2 optical counterpart at a redshift of 1.146 (ICowie et al.l 
2004). The MERLIN- VLA survey detected this radio source, 
finding a compact compone nt with a ~ 0.6 ar csec one-sided 
extension to the south-west (Mux low et al.l l2005). The infrared- 
bright object was also detected in the deep Chandra survey, in 
both the hard and soft bands. Our VLBI detects a compact core 
containing about half of the VLA flux density. 
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Fig. 10. VLA J123714+620823 



south-west (Muxlo w et al.l2005h : the radio structure is hosted by 
an / = 19.8 galaxy at z = 0.561 (Cowie 2012, priv.comm.). The 
VLBI observations detect a ~ 0.1 mJy compact r adio core in this 
(724 = -0.06 object. Recent VLA observations (Morriso n et al.l 
2010) indicate radio variability. 
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Fig. 12. VLA J123716+621733 
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3.2.8. VLA J1 23721 +621 129 

VLA J123721+621129 has no optical counterpart to / = 25 
and shows an inverted radio spectrum (a = -0.28) with a 
dominant compact radio com ponent and one-sided emission to 
the north ( Muxlow et al. 200 5|). Us i ng dee p imaging data from 
Barg er et all (1 1999b . iRichards et al.l (1 1999b report that this radio 
source is associated with a very red object having I - K > 5.2. 
Barg er et al.l (120081) list z — 1 .604. The compact radio compo- 
nent is slightly extended at ~ 25 mas resolution, and must be 
responsible for most if not all of the radio emission in this mod- 
erate radio-excess (^24 = -0.26) object. 



12 37 16.390 16.385 16.380 16.375 16.370 16.365 16.360 

RIGHT ASCENSION (J2000) 
Grey scale flux range- 38.40 86.00 MicroJY/BEAM 
Cont peak flux = 8.4124E-05 JY/BEAM 
Levs = 1.280E-05 " (-5, -3, 3, 5, 7, 9, 11, 13, 
15) 



Fig. 11. VLA J123716+621512 



4. Discussion 

This research has once more demonstrated the unique strength 
of high resolution radio astronomy in assessing the properties, 
i.e., mechanisms-at-work in the distant galaxy population. 
Our VLBI study bui lds upon the VLA+MERLIN studies of 
Muxl ow et al.l d2005l) . confirming, updating, and refining the 
results of that study. Being sensitive to compact radio emission 
having brightness temperatures in excess of 10 5 K, the VLBI 
observations unambiguously nail down AGN-related emission. 
Combining radio-spectral data with radio morphological and 
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i r 



(FR2-type) at z ~ 1-2 have radio luminosities in excess of 10 28 
W/Hz. The 5mJy FRl-type radio galaxy J123725+621 128 is 
identified with an / = 22.9 compact galaxy at unknown redshift; 
assuming z ~ 1 its radio luminosity would be ~ 10 25 W/Hz. 



62 1 1 30.05 — 



< 29.95 — 



± 



_L 



_L 



12 37 21.270 21.265 21.260 21.255 21.250 21.245 

RIGHT ASCENSION (J2000) 
Grey scale flux range- 46.8 105.0 MicroJY BEAM 
Cont peak flux = 1.0377E-04 JY/BEAM 
Levs = 1.560E-05 " (-5, -3, 3, 5, 7, 9, 11, 13, 
15) 



Fig. 13. VLA J123721+621129 



15ju infrared information, iMuxlow et al.l ([2005) classified 18 
out of their sample of 92 faint radio sources as AGN or AGN- 
candidate (AGNC). Their starburst and starburst-candidate class 
contains 48 objects. 28% of the sample, that is 26 objects, 
remain unclassified. It is furthermore very interesting to note 
that more than half of the total n=92 sample have Chandra 
X-ray counterparts, occurring both in the starburst and in the 
AGN class. With reference to the subfinal column of Table 1, 
our VLBI merely detects objects from the VLA+MERLIN 
AGN and AGNC classes, and not a single starburst object. 
Table 2 summarizes t he VLBI updat e with respect to the AGN 
and AGNC classes of IMuxlow et alj ((2005), and adds the new 
AGN. 

Of the 18 AGN and candidate AGN we confirm six 
AGN (J123608+621035, J123644+621133, J123646+621404, 
J123652+621444, J123714+620823, J123721 +621 1295), we 
nail down the AGN in four AGNC including the submm 
source GN16 (J123623+621642, J123700+620909/GN16, 
J123716+621512, J123716+621733) and in the z = 4.424 
starburst system with embedded AGN J123642+621331. In ad- 
dition we demonstrate the presence of an AGN in the z = 0.857 
object J123642+621545. It is seen that the VLBI misses out on 
seven faint (< 150/dy) AGN/ AGNC and on the relatively bright 
(5mJy) wide-angle-tail radio galaxy J123725+621 128, for the 
obvious reasons of sensitivity combined with resolution of ex- 
tended low surface brightness emission. Relevant (total) 1 .4GHz 
radio luminosities of the 1 1 VLBI detected radio sources with 
known redshifts, computed adopting a spectral shape S v oc v~° 5 
in the K-correction, range from ~ 1 x 10 23 to ~ 1 x 10 25 W/Hz, 
with outlier 2 x 10 26 W/Hz for the starburst with embedded 
AGN system J 126242+62 1331 at z = 4.424. This is in the 
range of radio-luminous Seyfert galaxies and edge-darkened 
double-lobed (FRl-type) radio galaxies (NGC 1068 (3C71) 
and M87 (Virgo A/3C 274) have P^ghz ~ 1 x 10 23 W/Hz) 
to radio-quiet QSOs. Note that large, double-lobed galaxies 



Considering their infrared/radio ratio 1724, the VLBI detected 
objects have at most a mild radio excess: the AGN generated 
radio emission is of similar strength as the starburst generated 
radio emission. Figure 14 illustrates that point. The diagram 
compares the (724-values of the 11 VLBI sources with known 
redshifts with the (724- vahies of the 2 4 starb urst objects with 
known redshift from IMuxlow et al.1 (120051) . Both samples 
contain upper limits, due to Spitzer upper limits. The diagram 
also specifies average (724-values for the samples (this time in- 
cluding the unknown redshift objects). While there is substantial 
overlap, the AGN stand out by a moderate (724-offset. The AGN 
sample displays an average (^24) = -0.12, whereas (^24) for 
SBs is 0.63+0. 43, which is consis tent with the previous result 
of 0.52+0.37 ( Beswick et all l2008h . Objects hav ing q 2 4 < 
are g enerally known as radio-excess objects (e.g. iDonlev et al.l 
2005). It is seen that four VLBI AGN, having q 2A = 0.08, 0.55, 
0.76 and 1.03, definitely would not classify as such. Others 
are seen to have a mild excess at most. Starburst candidate 
J123646+621629 at z — 0.502 is seen to have a radio excess 
(<724 = -0.7). It is puzzling that the VLBI observations do not 
detect this object; its MIPS upper limit may be erroneous or its 
radio core may be of variable strength (the object is unique in 
that sense). We stress that the integrated (VLA) 1.4 GHz flux 
densities often exceed the VLBI flux densities. The resolved, 
extended radio emission could originate in compact jets or 
lobes, or draw from diffuse galactic synchrotron radiation in 
the AGN host galaxies. The substantial 24yu emission detected 
in these objects would argue for an ISM origin in dusty host 
galaxies. The SED decomposition technique involving the NIR, 
M IR, FIR, submm ; and radio wavebands, recently employed 
bv lDel Moro et al.l (120121) . suggests symbiosis of dust-obscured 
star-formation and AGN in several dozens of weak, X-ray se- 
lected objects in the GOODS-N field. Their radio-excess sample 
contains eight of the VLBI AGN which provides confidence 
in both techniques. Surprisingly, a somewhat lower sSFR was 
established for these symbiotic objects, in comparison with 24/i 
select ed pure-starbur s ts. VL BI-undetected (candidate-) AGN 
from iDel Mor o et al ] (12012b have larger (724-values, indicating 
a smaller AGN contribution. Their radio excess flux densities 
are in the range of ~ 50 - 150^/Jy. The required global VLBI 
sensitivity (over a wide field) to tackle such radio sources is 
within reach; we are currently planning such observations. A 
detailed comparison of the results of the SED decomposition 
and the radio imaging techniques is beyond the scope of the 
present analysis and is postponed to a later paper. 

Whereas X-ray selection is a powerful tool to select 
AGN samples, the present research underlines its shortcom- 
ings. Two of the VLBI-selected AGN, J 123623 +62 1642 and 
J123700+620909/GN16 remai n without X-ray detect ion. The 
2Ms Chandra observations of lAlexander et all (120031) yielded 
upper limits logL2-iokeV < 42.6 erg/s for both objects. The 
dusty nature of their host galaxies, referred to above, could be 
responsible for X-ray attenuation; this issue again wilt be dealt 
with in more detail in a later paper (the X-ray luminosities for 
the detected VLBI AGN are in the range 10 41 - 10 44 erg/s). It 
is in the meantime clear that combination of arcsec scale and 
milliarcsec scale radio imaging provides a very powerful and 
in principle straightforward tool to select samples of AGN and 
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star-forming galaxies among the distant galaxy population. It is 
no surprise that such objects are among the main science drivers 
of the LOFAR and SKA radio telescopes. 

To summarize, faint AGN are found in a variety of star- 
forming galaxies, with a wide range of dust obscurations, and 
a relative occurrence of ~ 25%. A subsequent paper will address 
the properties as well as selection techniques of faint AGN in 
more detail, and discuss the general AGN evolution over cosmic 
time. 



5. Conclusions 

Deep, high-resolution global VLBI observations at milliarcsec 
resolution together with deep arcsec-scale resolution radio imag- 
ing provide an extremely powerful way to separate faint AGN 
and starburst galaxies, and assess their symbiotic occurrence 
among the distant galaxy population. Using these techniques on 
the Hubble Deep Field North and its flanking fields, it is found 
that about 25% of star-forming galaxies of various types harbour 
faint AGN. Some classify as starbursts, such as the z - 4.42 sys- 
tem J123642+621331 and the z = 1.68 submm galaxy GN16/ 
J 123700+620909, and some can be X-ray obscured. 
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rap 





Fig. 1. Composite image of the radio (WSRT 1 .4 GHz)-optical overlay image of the HDF-N and HFF, surrounded by postage stamp 
images of the twelve compact VLBI-detected radio sources. The cyan circles represent annuli of decreasing resolution and sensitivity, 
and are drawn at 2, 4, 6, and 8 arcmin radius w.r.t. the phase center which coincides with radio AGN VLA J123642+621331 (see 
text Section 3). 
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Fig. 14. Values of ^24 for SBs (black empty squares) and 1 1 global VLBI detections (red filled squares), as function of redshift. The thick dashed 
(black) and solid (red) lines represent the mean values for each sample, respectively, and the Icr deviation is shown by thinner black dashed and 
red lines. 
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Table 1. Summary: compact radio sources in HDF-N and HFF. 



Name 


R.A. 


Dec. 


z 


S VIA " 


S WSRT " 


S Global VLB1 


r b 




Radio 


Classification 


924 




[hh mm ss] 


[+dd mm ss] 




b"Jy] 




L"Jy] 




[lO^WirT 2 ] 


structure 






Inner HDF-N and adjacent HFF 


J123642+621331 (core) 
J123642+621331 (jet) 


12 36 42.0908 


62 13 31.425 


4.424 e 


467 


489/489 


227+26 
121+21 


209±5.59 


0.28 


C1E 


SB+AGN 


-0.35 


J123644+621133 


12 36 44.3870 


62 11 33.145 


1.050 f 


1290 


1190/1606 


309+27 


< 80 


0.24 


FR1 


AGN (FR1) 


< -1.21 


J123646+621404 


12 36 46.3321 


62 14 04.693 


0.961 s 


179 


187/187 


247+14 


213±5.22 


24.6 


C1E 


AGN 


0.08 


J123652+621444 


12 36 52.8839 


62 14 44.076 


0.321 * 


168 


237/276 


83±13 


< 80 


0.46 


C1E 


AGN 


< -0.32 


Outer HFF 


J123608+621035 


12 36 08.1195 


62 10 35.898 


0.681 s 


217 


190/246 


140+30 


2300±13.50 


1.89 


CE 


AGN 


1.03 


J123623+621642 


12 36 23.5436 


62 16 42.754 


1.918 '' 


481 


476/476 


327+50 


< 80 




C1E 


AGN candidate 


< -0.78 


J123642+621545 


12 36 42.2123 


62 15 45.521 


0.857 ' 


150 


88/89 


343±101 


866±10.30 


2.46 


CE 


U (SB+AGN) 


0.76 


J123700+620909 


12 37 00.2480 


62 09 09.778 


1.68 ^ 


324 


236/236 


147±34 


267±7.62 




C1E 


AGN candidate (GN16) 


-0.08 


J123714+620823 


12 37 14.9414 


62 08 23.208 


(0.847) * 


1350 


1853/1853 


645±80 


215±6.67 


0.6 


C 


AGN 


-0.80 


J123716+621512 


12 37 16.3740 


62 15 12.343 


0.561 ' 


187 


166/166 


110+25 


162±9.40 


0.23 


C1E 


AGN candi. (AGN) 


-0.06 


J123716+621733 


12 37 16.6811 


62 17 33.327 


1.146 s 


346 


362/362 


177+25 


1240± 15.50 


21.69 


C1E 


AGN candi. (AGN) 


0.55 


J123721+621129 


12 37 21.2539 


62 11 29.954 


1.56'" 


383 


381/381 


254+51 


(210+6.29) d 


0.51 


C1E 


AGN 


-0.26 



" iMuxlow et alj j2005h: the WSRT entries represent peak/total values • 

b Spitzer. merged IRAC a nd MIPS GOODS-N catalog (R. Chary - 201 1 private communication) 9 

c lAlexander et all d2003h : lDonlev et al.l d2005h g 

d ~ 1.5" ofFsetfrom MIPS position, ~ 1" offset from IRAC & 

e IWaddington et al.Ul999l ) g 

^ ICohen et al.U2000h D 

s Cowie et al.H2004l) ^ 

'' IChapman etai] d2004l) g 

' Hornschemeier et al. d200lh o- 

J lPopeet alJ J2006b | 

* L. Cowie - 2010 private communication: tentative value ^ 
1 L. Cowie - 2012 pri vate communication; not to be confused with J 1237 16+62 151 1 at z = 0.231 
m lBargeretaHd2002l) 
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Table 2. VLBI detected and undetected AGN and AGN candidates from the VLA+MERLIN 1.4 GHz survey dMuxlow et all2 005). 
Note that the VLBI in addition detected the starburst+AGN system J123642+621331 at z = 4.42 and the unclear identification 
J123642+621545 at z = 0.857. The VLBI sensitivity ranges from ~ 50yuJy in the field center to ~ 150yuJy at r p ha Sece nter > 6 arcmin 
(see text). 



Cat.number" 


Name 6 


Alt.name 


1.4 GHz VLA flux density 


VLA class c 


VLBI detection 


4 


J123608+621035 




217/Jy 


AGN 


yes 


15 


J123618+621635 




47//Jy 


AGNC 


no 


18 


J123620+620844 




123/iJy 


AGN 


no 


21 


J123622+621544 




84pJy 


AGNC 


no 


24 


J123623+621642 




481/Jy 


AGNC 


yes 


39 


J123640+621009 




87/Jy 


AGNC 


no 


43 


J123644+621133 




1290^ 


AGN 


yes 


47 


J123646+621404 




179//Jy 


AGN 


yes 


56 


J123652+621444 




168/Jy 


AGN 


yes 


64 


J123700+620909 


GN16 


324/Jy 


AGNC 


yes 


72 


J123709+620837 




72pJy 


AGN 


no 


73 


J123709+620841 




68//Jy 


AGNC 


no 


74 


J123711+621330 




132/Jy 


AGNC 


no 


78 


J123714+620823 




1350yuJy 


AGN 


yes 


79 


J123716+621512 




187^/Jy 


AGNC 


yes 


81 


J123716+621733 




346/iJy 


AGNC 


yes 


85 


J123721+621129 




383/iJy 


AGN 


yes 


89 


J123725+621128 




5mJy 


AGN 


no 



" number in VLA/MERLIN catalog I Muxlow et al. 2005), at increasing RA 
* J2 000 name of radio sou rce dMuxlow et al.ll2005h 

c the lMuxlow et a l. (2005) classification combines radio and infrared (ISO) information 
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